Rigorous Screened Interactions for Realistic Correlated Electron Systems.

We derive a widely applicable first-principles approach for determining two-body, static effective interactions for low-energy Hamiltonians with quantitative accuracy. The algebraic construction rigorously conserves all instantaneous two-point correlation functions in a chosen model space at the level of the random phase approximation, improving upon the traditional uncontrolled static approximations. Applied to screened interactions within a quantum embedding framework, we demonstrate these faithfully describe the relaxation of local subspaces via downfolding high-energy physics in molecular systems, as well as enabling a systematically improvable description of the long-range plasmonic contributions in extended graphene.

Effective Reconstruction of Expectation Values from Ab Initio Quantum Embedding.

Quantum embedding is an appealing route to fragment a large interacting quantum system into several smaller auxiliary "cluster" problems to exploit the locality of the correlated physics. In this work, we critically review approaches to recombine these fragmented solutions in order to compute nonlocal expectation values, including the total energy. Starting from the democratic partitioning of expectation values used in density matrix embedding theory, we motivate and develop a number of alternative approaches, numerically demonstrating their efficiency and improved accuracy as a function of increasing cluster size for both energetics and nonlocal two-body observables in molecular and solid state systems. These approaches consider the N-representability of the resulting expectation values via an implicit global wave function across the clusters, as well as the importance of including contributions to expectation values spanning multiple fragments simultaneously, thereby alleviating the fundamental locality approximation of the embedding. We clearly demonstrate the value of these introduced functionals for reliable extraction of observables and robust and systematic convergence as the cluster size increases, allowing for significantly smaller clusters to be used for a desired accuracy compared to traditional approaches in ab initio wave function quantum embedding.

A "moment-conserving" reformulation of GW theory.

We show how to construct an effective Hamiltonian whose dimension scales linearly with system size, and whose eigenvalues systematically approximate the excitation energies of GW theory. This is achieved by rigorously expanding the self-energy in order to exactly conserve a desired number of frequency-independent moments of the self-energy dynamics. Recasting GW in this way admits a low-scaling O[N4] approach to build and solve this Hamiltonian, with a proposal to reduce this further to O[N3]. This relies on exposing a novel recursive framework for the density response moments of the random phase approximation, where the efficient calculation of its starting point mirrors the low-scaling approaches to compute RPA correlation energies. The frequency integration of GW, which distinguishes so many different GW variants, can be performed without approximation directly in this moment representation. Furthermore, the solution to the Dyson equation can be performed exactly, avoiding analytic continuation, diagonal approximations, or iterative solutions to the quasiparticle equation, with the full-frequency spectrum obtained from the complete solution of this effective static Hamiltonian. We show how this approach converges rapidly with respect to the order of the conserved self-energy moments and is applied across the GW100 benchmark dataset to obtain accurate GW spectra in comparison to traditional implementations. We also show the ability to systematically converge all-electron full-frequency spectra and high-energy features beyond frontier excitations, as well as avoiding discontinuities in the spectrum, which afflict many other GW approaches.

Constructing "Full-Frequency" Spectra via Moment Constraints for Coupled Cluster Green's Functions.

We propose an approach to build "full-frequency" quasiparticle spectra from conservation of a set of static expectation values. These expectation values define the moments of the spectral distribution, resulting in an efficient and systematically improvable expansion. By computing these initial moment constraints at the coupled-cluster level, we demonstrate convergence in both correlated state-specific and full spectral quantities, while requiring a fraction of the effort of traditional Green's function approaches. Tested across the GW100 benchmark set for charged excitation spectra, we can converge frontier excitations to within the inherent accuracy of the CCSD approximation, while obtaining a simultaneous representation of the entire excitation spectrum at all energy scales.

Scalable and Predictive Spectra of Correlated Molecules with Moment Truncated Iterated Perturbation Theory.

A reliable and efficient computation of the entire single-particle spectrum of correlated molecules is an outstanding challenge in the field of quantum chemistry, with standard density functional theory approaches often giving an inadequate description of excitation energies and gaps. In this work, we expand upon a recently introduced approach that relies on a fully self-consistent many-body perturbation theory coupled to a nonperturbative truncation of the effective dynamics at each step. We show that this yields a low-scaling and accurate method across a diverse benchmark test set that is capable of treating moderate levels of strong correlation effects, and we detail an efficient implementation for applications involving up to ∼1000 orbitals on parallel resources. We then use this method to characterize the spectral properties of the antimalarial drug molecule artemisinin, resolving discrepancies in previous works concerning the active sites of the lowest-energy fundamental excitations of the system.

Equation of state of atomic solid hydrogen by stochastic many-body wave function methods.

We report a numerical study of the equation of state of crystalline body-centered-cubic (BCC) hydrogen, tackled with a variety of complementary many-body wave function methods. These include continuum stochastic techniques of fixed-node diffusion and variational quantum Monte Carlo and the Hilbert space stochastic method of full configuration-interaction quantum Monte Carlo. In addition, periodic coupled-cluster methods were also employed. Each of these methods is underpinned with different strengths and approximations, but their combination in order to perform reliable extrapolation to complete basis set and supercell size limits gives confidence in the final results. The methods were found to be in good agreement for equilibrium cell volumes for the system in the BCC phase.

Four-component full configuration interaction quantum Monte Carlo for relativistic correlated electron problems.

An adaptation of the full configuration interaction quantum Monte Carlo (FCIQMC) method is presented for correlated electron problems containing heavy elements and the presence of significant relativistic effects. The modified algorithm allows for the sampling of the four-component spinors of the Dirac-Coulomb(-Breit) Hamiltonian within the relativistic no-pair approximation. The loss of spin symmetry and the general requirement for complex-valued Hamiltonian matrix elements are the most immediate considerations in expanding the scope of FCIQMC into the relativistic domain, and the alternatives for their efficient implementation are motivated and demonstrated. For the canonical correlated four-component chemical benchmark application of thallium hydride, we show that the necessary modifications do not particularly adversely affect the convergence of the systematic (initiator) error to the exact correlation energy for FCIQMC calculations, which is primarily dictated by the sparsity of the wavefunction, allowing the computational effort to somewhat bypass the formal increases in Hilbert space dimension for these problems. We apply the method to the larger problem of the spectroscopic constants of tin oxide, correlating 28 electrons in 122 Kramers-paired spinors, finding good agreement with experimental and prior theoretical relativistic studies.

A Bayesian inference framework for compression and prediction of quantum states.

The recently introduced Gaussian Process State (GPS) provides a highly flexible, compact, and physically insightful representation of quantum many-body states based on ideas from the zoo of machine learning approaches. In this work, we give a comprehensive description of how such a state can be learned from given samples of a potentially unknown target state and show how regression approaches based on Bayesian inference can be used to compress a target state into a highly compact and accurate GPS representation. By application of a type II maximum likelihood method based on relevance vector machines, we are able to extract many-body configurations from the underlying Hilbert space, which are particularly relevant for the description of the target state, as support points to define the GPS. Together with an introduced optimization scheme for the hyperparameters of the model characterizing the weighting of modeled correlation features, this makes it possible to easily extract physical characteristics of the state such as the relative importance of particular correlation properties. We apply the Bayesian learning scheme to the problem of modeling ground states of small Fermi-Hubbard chains and show that the found solutions represent a systematically improvable trade-off between sparsity and accuracy of the model. Moreover, we show how the learned hyperparameters and the extracted relevant configurations, characterizing the correlation of the wave function, depend on the interaction strength of the Hubbard model and the target accuracy of the representation.

Efficient Excitations and Spectra within a Perturbative Renormalization Approach.

We present a self-consistent approach for computing the correlated quasiparticle spectrum of charged excitations in iterative O[N5] computational time. This is based on the auxiliary second-order Green's function approach [Backhouse, O. J. Chem. Theory Comput., 2000], in which a self-consistent effective Hamiltonian is constructed by systematically renormalizing the dynamical effects of the self-energy at second-order perturbation theory. From extensive benchmarking across the W4-11 molecular test set, we show that the iterative renormalization and truncation of the effective dynamical resolution arising from the 2h1p and 1h2p spaces can substantially improve the quality of the resulting ionization potential and electron affinity predictions compared to benchmark values. The resulting method is shown to be superior in accuracy to similarly scaling quantum chemical methods for charged excitations in EOM-CC2 and ADC(2), across this test set, while the self-consistency also removes the dependence on the underlying mean-field reference. The approach also allows for single-shot computation of the entire quasiparticle spectrum, which is applied to the benzoquinone molecule and demonstrates the reduction in the single-particle gap due to the correlated physics, and gives direct access to the localization of the Dyson orbitals.

Recent developments in the PySCF program package.

PySCF is a Python-based general-purpose electronic structure platform that supports first-principles simulations of molecules and solids as well as accelerates the development of new methodology and complex computational workflows. This paper explains the design and philosophy behind PySCF that enables it to meet these twin objectives. With several case studies, we show how users can easily implement their own methods using PySCF as a development environment. We then summarize the capabilities of PySCF for molecular and solid-state simulations. Finally, we describe the growing ecosystem of projects that use PySCF across the domains of quantum chemistry, materials science, machine learning, and quantum information science.

NECI: N-Electron Configuration Interaction with an emphasis on state-of-the-art stochastic methods.

We present NECI, a state-of-the-art implementation of the Full Configuration Interaction Quantum Monte Carlo (FCIQMC) algorithm, a method based on a stochastic application of the Hamiltonian matrix on a sparse sampling of the wave function. The program utilizes a very powerful parallelization and scales efficiently to more than 24 000 central processing unit cores. In this paper, we describe the core functionalities of NECI and its recent developments. This includes the capabilities to calculate ground and excited state energies, properties via the one- and two-body reduced density matrices, as well as spectral and Green's functions for ab initio and model systems. A number of enhancements of the bare FCIQMC algorithm are available within NECI, allowing us to use a partially deterministic formulation of the algorithm, working in a spin-adapted basis or supporting transcorrelated Hamiltonians. NECI supports the FCIDUMP file format for integrals, supplying a convenient interface to numerous quantum chemistry programs, and it is licensed under GPL-3.0.

Driven Imposters: Controlling Expectations in Many-Body Systems.

We present a framework for controlling the observables of a general correlated electron system driven by an incident laser field. The approach provides a prescription for the driving required to generate an arbitrary predetermined evolution for the expectation value of a chosen observable, together with a constraint on the maximum size of this expectation. To demonstrate this, we determine the laser fields required to exactly control the current in a Fermi-Hubbard system under a range of model parameters, fully controlling the nonlinear high-harmonic generation and optically observed electron dynamics in the system. This is achieved for both the uncorrelated metalliclike state and deep in the strongly correlated Mott insulating regime, flipping the optical responses of the two systems so as to mimic the other, creating "driven imposters." We also present a general framework for the control of other dynamical variables, opening a new route for the design of driven materials with customized properties.

Efficient and stochastic multireference perturbation theory for large active spaces within a full configuration interaction quantum Monte Carlo framework.

Full Configuration Interaction Quantum Monte Carlo (FCIQMC) has been effectively applied to very large configuration interaction (CI) problems and was recently adapted for use as an active space solver and combined with orbital optimization. In this work, we detail an approach within FCIQMC to allow for efficient sampling of fully internally contracted multireference perturbation theories within the same stochastic framework. Schemes are described to allow for the close control over the resolution of stochastic sampling of the effective higher-body intermediates within the active space. It is found that while complete active space second-order perturbation theory seems less amenable to a stochastic reformulation, strongly contracted N-Electron Valence second-order Perturbation Theory (NEVPT2) is far more stable, requiring a similar number of walkers to converge the sc-NEVPT2 expectation values as to converge the underlying CI problem. We demonstrate the application of the stochastic approach to the computation of sc-NEVPT2 within a (24, 24) active space in a biologically relevant system and show that small numbers of walkers are sufficient for a faithful sampling of the sc-NEVPT2 energy to chemical accuracy, despite the active space already exceeding the limits of practicality for traditional approaches. This raises prospects of an efficient stochastic solver for multireference chemical problems requiring large active spaces, with an accurate treatment of external orbitals.

Wave Function Perspective and Efficient Truncation of Renormalized Second-Order Perturbation Theory.

We present an approach to a renormalized second-order Green's function perturbation theory (GF2), which avoids all dependency on continuous variables, grids, or explicit Green's functions and is instead formulated entirely in terms of static quantities and wave functions. Correlation effects from MP2 diagrams are iteratively incorporated to modify the underlying spectrum of excitations by coupling the physical system to fictitious auxiliary degrees of freedom, allowing for single-particle orbitals to delocalize into this additional space. The overall approach is shown to be rigorously O[N5], after an appropriate compression of this auxiliary space. This is achieved via a novel scheme, which ensures that a desired number of moments of the underlying occupied and virtual spectra are conserved in the compression, allowing a rapid and systematically improvable convergence to the limit of the effective dynamical resolution. The approach is found to then allow for the qualitative description of stronger correlation effects, avoiding the divergences of MP2, as well as its orbital-optimized version. On application to the G1 test set, we find that modification up to only the third spectral moment of the underlying spectrum from which the double excitations are built are required for accurate energetics, even in strongly correlated regimes. This is beyond simple self-consistent changes to the density matrix of the system but far from requiring a description of the full dynamics of the frequency-dependent self-energy.

Energy-weighted density matrix embedding of open correlated chemical fragments.

We present a multiscale approach to efficiently embed an ab initio correlated chemical fragment described by its energy-weighted density matrices and entangled with a wider mean-field many-electron system. This approach, first presented by Fertitta and Booth [Phys. Rev. B 98, 235132 (2018)], is here extended to account for realistic long-range interactions and broken symmetry states. The scheme allows for a systematically improvable description in the range of correlated fluctuations out of the fragment into the system, via a self-consistent optimization of a coupled auxiliary mean-field system. It is discussed that the method has rigorous limits equivalent to the existing quantum embedding approaches of both dynamical mean-field theory and density matrix embedding theory, to which this method is compared, and the importance of these correlated fluctuations is demonstrated. We derive a self-consistent local energy functional within the scheme and demonstrate the approach for hydrogen rings, where quantitative accuracy is achieved despite only a single atom being explicitly treated.

Response Formalism within Full Configuration Interaction Quantum Monte Carlo: Static Properties and Electrical Response.

We formulate a general, arbitrary-order stochastic response formalism within the full configuration interaction quantum Monte Carlo framework. This modified stochastic dynamic allows for the exact response properties of correlated multireference electronic systems to be systematically converged upon for systems far out of reach of traditional exact treatments. This requires a simultaneous coupled evolution of a response state alongside the zero-order state, which is shown to be stable, nontransient, and unbiased. We demonstrate this with application to the static dipole polarizability of molecular systems and, in doing so, resolve a discrepancy between restricted and unrestricted high-level coupled-cluster linear response results which were the high-accuracy benchmark in the literature.

Density matrices in full configuration interaction quantum Monte Carlo: Excited states, transition dipole moments, and parallel distribution.

We present developments in the calculation of reduced density matrices (RDMs) in the full configuration interaction quantum Monte Carlo (FCIQMC) method. An efficient scheme is described to allow storage of RDMs across distributed memory, thereby allowing their calculation and storage in large basis sets. We demonstrate the calculation of RDMs for general states by using the recently introduced excited-state FCIQMC approach [N. S. Blunt et al., J. Chem. Phys. 143, 134117 (2015)] and further introduce calculation of transition density matrices in the method. These approaches are combined to calculate excited-state dipole and transition dipole moments for heteronuclear diatomic molecules, including LiH, BH, and MgO, and initiator error is investigated in these quantities.

A comparison between quantum chemistry and quantum Monte Carlo techniques for the adsorption of water on the (001) LiH surface.

We present a comprehensive benchmark study of the adsorption energy of a single water molecule on the (001) LiH surface using periodic coupled cluster and quantum Monte Carlo theories. We benchmark and compare different implementations of quantum chemical wave function based theories in order to verify the reliability of the predicted adsorption energies and the employed approximations. Furthermore we compare the predicted adsorption energies to those obtained employing widely used van der Waals density-functionals. Our findings show that quantum chemical approaches are becoming a robust and reliable tool for condensed phase electronic structure calculations, providing an additional tool that can also help in potentially improving currently available van der Waals density-functionals.

Projector Quantum Monte Carlo Method for Nonlinear Wave Functions.

We reformulate the projected imaginary-time evolution of the full configuration interaction quantum Monte Carlo method in terms of a Lagrangian minimization. This naturally leads to the admission of polynomial complex wave function parametrizations, circumventing the exponential scaling of the approach. While previously these functions have traditionally inhabited the domain of variational Monte Carlo approaches, we consider recent developments for the identification of deep-learning neural networks to optimize this Lagrangian, which can be written as a modification of the propagator for the wave function dynamics. We demonstrate this approach with a form of tensor network state, and use it to find solutions to the strongly correlated Hubbard model, as well as its application to a fully periodic ab initio graphene sheet. The number of variables which can be simultaneously optimized greatly exceeds alternative formulations of variational Monte Carlo methods, allowing for systematic improvability of the wave function flexibility towards exactness for a number of different forms, while blurring the line between traditional variational and projector quantum Monte Carlo approaches.

From plane waves to local Gaussians for the simulation of correlated periodic systems.

We present a simple, robust, and black-box approach to the implementation and use of local, periodic, atom-centered Gaussian basis functions within a plane wave code, in a computationally efficient manner. The procedure outlined is based on the representation of the Gaussians within a finite bandwidth by their underlying plane wave coefficients. The core region is handled within the projected augment wave framework, by pseudizing the Gaussian functions within a cutoff radius around each nucleus, smoothing the functions so that they are faithfully represented by a plane wave basis with only moderate kinetic energy cutoff. To mitigate the effects of the basis set superposition error and incompleteness at the mean-field level introduced by the Gaussian basis, we also propose a hybrid approach, whereby the complete occupied space is first converged within a large plane wave basis, and the Gaussian basis used to construct a complementary virtual space for the application of correlated methods. We demonstrate that these pseudized Gaussians yield compact and systematically improvable spaces with an accuracy comparable to their non-pseudized Gaussian counterparts. A key advantage of the described method is its ability to efficiently capture and describe electronic correlation effects of weakly bound and low-dimensional systems, where plane waves are not sufficiently compact or able to be truncated without unphysical artifacts. We investigate the accuracy of the pseudized Gaussians for the water dimer interaction, neon solid, and water adsorption on a LiH surface, at the level of second-order Møller-Plesset perturbation theory.